Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 30 March 2010 (MN KTeX style file v2.2) 



Feedback under the microscope II: heating, gas uphft, and mixing in 
the nearest cluster core 

N. Werner^*, A. Simionescu^t, E. T. Million\ S. W. Allen\ P. E. J. Nulsen^, 

A. von der Linden^ S. M. Hansen^, H. Bohringer'^, E. Churazov^'^, A. C. Fabian^, 

W. R. Forman^, C. Jones^, J. S. Sanders^, and G. B. Taylor^ 

^Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305-4060, USA 
^ and SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA 
^Harvard- Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 

^University of California Observatories & Department of Astronomy, University of California, Santa Cruz, CA 95064, USA 
^ ^ Max-Planck-Institut fUr extraterrestrische Physik, Giessenbachstr, 85748 Garching, Germany 
^ Max-Planck-Institut fUr Astrophysik, Karl-Schwarzschild-Strasse 1, 85741 Garching, Germany 
^ Space Research Institute (IKI), Profsoyznaya 84/32, Moscow 117810, Russia 
^Institute of Astronomy, Madingley Road, Cambridge CB3 OHA 
QQ ^Department of Physics and Astronomy, University of New Mexico, Albuquerque, NM 87131, USA 

O 30 March 2010 

u 

ABSTRACT 

Using a combination of deep (574 ks) Chandra data, XMM-Newton high-resolution spectra, 
Q and optical Ha+[Nii] images, we study the nature and spatial distribution of the multiphase 

^ plasma in M 87. Our results provide direct observational evidence of 'radio mode' AGN feed- 

C/^ back in action, stripping the central galaxy of its lowest entropy gas and therefore preventing 

star-formation. This low entropy gas was entrained with and uplifted by the buoyantly ris- 
ing relativistic plasma, forming long "arms". A number of arguments suggest that these arms 
are oriented within 15°-30° of our line-of-sight. The mass of the uplifted gas in the arms is 
^ comparable to the gas mass in the approximately spherically symmetric 3.8 kpc core, demon- 

strating that the AGN has a profound effect on its immediate surroundings. The coolest X-ray 
emitting gas in M 87 has a temperature of -0.5 keV and is spatially coincident with Ha-h[N ii] 
CO nebulae, forming a multiphase medium where the cooler gas phases are arranged in magne- 

tized filaments. We place strong upper limits of 0.06 Mo/yr (at 95 per cent confidence) on the 
amount of plasma cooling radiatively from 0.5 to 0.25 keV and show that a uniform, volume- 
averaged heating mechanism could not be preventing the cool gas from further cooling. All 
of the bright Ha filaments in M 87 appear in the downstream region of the < 3 Myr old 
i-H shock front, at smaller radii than -0.6'. We suggest that shocks induce shearing around the 

^ filaments, thereby promoting mixing of the cold gas with the ambient hot ICM via instabil- 

• i-H ities. By bringing hot thermal particles into contact with the cool, line-emitting gas, mixing 

can supply the power and ionizing particles needed to explain the observed optical spectra. 
^ Furthermore, mixing of the coolest X-ray emitting plasma with the cold optical line emitting 

filamentary gas promotes efficient conduction between the two phases, allowing non-radiative 
cooling which could explain the lack of X-ray gas with temperatures under 0.5 keV. 

Key words: X-rays: galaxies: clusters - galaxies: individual: M 87 - galaxies: intergalactic 
medium - cooling flows 



1 INTRODUCTION 

If the energy radiated away at the centres of so-called "cool-core" 
clusters of galaxies, which show sharp X-ray surface brightness 

* Chandra/Einstein fellow, E-mail: norbertw@stanford.edu 
t Einstein fellow 



peaks and central temperature dips, came only from the thermal en- 
ergy of the hot, diffuse intra-cluster medium (ICM), the ICM would 
cool and form stars at rates orders of magnitude above what the 
observations suggest (see Peterson & Fabianp OOG' for a review). 
It is currently believed that the energy which offsets the cooling 
is provided by the interaction between the active galactic nuclei 
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(AGN) in the central dominant galaxies and the ICM (see |Chu-| 
Irazov et aL|[2000l |200T] [20021 [McNamara & Nulsen|[2007l for a 
review). Through a tight feedback loop, it is thought that the AGN 
can provide enough energy to prevent catastrophic cooling. Recent 
deep observations of nearby bright cooling core clusters such as 
Perseus, M 87, Centaurus, and Hydra A revealed AGN induced 
weak shocks and sound waves with sufficient energy flux to coun- 
terbalance the radiative cooling ("Sanders & Fabian"2007', "Forman 
[eta l. 2005 , Sa nders & Fabian 2008 , Nulsen et al. 2005 , Simionescu 
|et a l. 2009a). To the first order it is thus understood where the en- 
ergy comes from and how it gets transported. However, the detailed 
physics of the AGN/ICM interaction is not yet clear. 

M 87, the central dominant galaxy of the Virgo cluster, at a 
distance of only 16.1 Mpc ( Tonry et al. 2001 ) is the nearest. X-ray 
bright cool core. It is an ideal system for detailed studies of the en- 
ergy input from the AGN to the hot, cooling ICM. After the Perseus 
cluster, M 87 is the second brightest extended extragalactic object 
in soft X-rays. The most striking X-ray features in M 87 are its two 
"X-ray arms", first detected by Feigelson et al. (1987 ) using Ein- 
stein Observatory data and later studied in detail by Bo hringer et al.| 
( ,1995| ) using ROSAT and VLA radio observations. The X-ray arms 
extend to the East and Southwest from the centre of the galaxy and 
they spatially correlate with the radio emission. [Bohringer et al.| 
( J^95 ) found that the X-ray arms are due to cooler gas, possibly 
uplifted from the centre of the galaxy. These results and the sub- 
sequent high-quality radio data by jOwen et aT] ( |2000| ) led [Chura-[ 
jzov et al.[ ( [2001| ) to argue that the X-ray and radio morphology can 
be explained by bubbles of radio-emitting plasma rising buoyantly 
through the hot ICM. These rising bubbles entrain and uplift cooler 
gas from the centre of the galaxy, which then further adiabatically 
cools. [Churazov et al.[ ( [200 1| ) suggested that the buoyantly rising 
bubbles dissipate energy into sound waves, internal waves, turbu- 
lent motion in the wake, and potential energy of the uplifted gas, 
all of which could provide heating to the cooling core region (see 
also e.g.jBriiggen & Kaiser|20 02, Kai ser|2003 [ [Bruggen|2003 1 [De[ 
|Young|2003||Ruszkowskietri.|2004a|b[[Heinz&Churazov|2005| ) 

Observations with Chandra and XMM-Newton greatly en- 
hanced our view of M 87. Using XMM-Newton observations, Bel- 
jsole et al.[ ( [200T] ), [Matsushita et al.j ( |2002l ), and .Molendi, ( ,2002j ) 
showed that the regions associated with radio arms have two- 
temperature (in the ranges of kT ~ 0.8-1 keV and kT ~ 1.6- 
2.5 keV) or multi-temperature structure. However, plasma below 
0.8 keV is largely absent ( [Bohringer et al.[[2002[ [SakelHou et al.[ 
|2002| ). [Young et al.[ ( |200 2') analyzed an early, relatively short (40 ks) 
Chandra observation which showed cavities and edges in the hot 
plasma. They compared the X-ray morphology with the 6 cm ra- 
dio ( [Hines et al. [19 89 ) and HQf+[Nii] emission (Sparks etal.|1993| ) 
and showed that the optical filaments lie outside of X-ray cavities 
filled with radio plasma; some are along the edges of the cavities, 
and several optical filaments coincide with knots of cooler X-ray 
gas. Using a longer (~150 ks) Chandra observation. Sparks et al. 
( [2004^ studied the relation of the X-ray and optical filaments and 
concluded that electron conduction from the hot X-ray emitting 
plasma to the cooler phase provides a quantitatively acceptable 
energy source for the optical filaments. Using the same Chandra 
data, [Forman et al.[ ( [2005| ) identified shock fronts associated with 
an AGN outburst about 1-2x10^ yr ago. They argued that shock 
fronts may be the most significant channel for heating (i.e., entropy 
increase) of the ICM near to the AGN. 

Subsequent Chandra observations increased the total expo- 
sure time to 500 ks, which allowed [Forman et aL] ( [2007[ ) to resolve 
a web of filamentary structure in both X-ray arms. Both arms show 



Table 1. Summary of the Chandra observations. Colums list the observation 
ID, detector, observation mode, exposure after cleaning, and observation 
date. 



Obs. ID 


Detector 


Mode 


Exposure (ks) 


Obs. date 


2707 


ACIS-S 


FAINT 


82.9 


Jul. 6 2002 


3717 


ACIS-S 


FAINT 


11.1 


Jul. 5 2002 


5826 


ACIS-I 


VFAINT 


126.8 


Mar. 3 2005 


5827 


ACIS-I 


VFAINT 


156.2 


May 5 2005 


5828 


ACIS-I 


VFAINT 


33.0 


Nov. 17 2005 


6186 


ACIS-I 


VFAINT 


51.5 


Jan. 31 2005 


7210 


ACIS-I 


VFAINT 


30.7 


Nov. 16 2005 


7211 


ACIS-I 


VFAINT 


16.6 


Nov. 16 2005 


7212 


ACIS-I 


VFAINT 


65.2 


Nov 14 2005 



narrow filaments, with a length-to-width ratio of up to ~50. How- 
ever, while the Southwestern arm has only a single set of filaments, 
the Eastern arm shows multiple sets of filaments and bubbles, in- 
creasing in scale with the distance from the centre. The high spatial 
resolution of Chandra also highlights the fact that while the East- 
ern X-ray arm is cospatial with the radio arm, the radio emission 
bends around and avoids the Southwestern X-ray arm. A residual 
image of M 87 showing the inner cavities and the X-ray arms ob 
tained using all Chandra ACIS-I data is shown in Fig.[T](from Mil 
lion et al.][2010| ), together with the 90 cm radio image by ,Owen 
et al.[pOOOj ). The temperature structure and metallicity of the X- 
ray arms was studied in detail by [Simionescu et al.[ ( [2008| ) using 
deep (120 ks) XMM-Newton data. These authors confirmed that the 
arms are multiphase with a temperature distribution between 0.6- 
3.2 keV. This temperature distribution is also consistent with that 
found by spectral fits to the XMM-Newton Reflection Grating Spec- 
trometer (RGS) data ( Werner et al.[2006|), which show the presence 
of weak Fe xvii lines. Simionescu et al.| (g008 ) estimate the total 
mass of gas below 1 .5 keV, most of which was uplifted by the AGN, 
as 5 X 10^ Mo. 

This is the second in a series of papers (following the work 
by [Million et al.|2 (}T0 hereafter Paper I) presenting detailed spa- 
tially resolved spectroscopy of M 87, exploiting the superb spa- 
tial resolution and excellent photon statistics of the deep (574 ks) 
Chandra observation, to study AGN feedback and ICM physics. In 
this paper, we study for the first time the detailed nature and the 
spatial distribution of the multiphase plasma associated with the 
X-ray arms. Chandra allows us to study detailed structure in the 
cooler phases at a much higher spatial resolution than is possible 
with XMM-Newton. Our data allow us to sample well the small 
substructure seen in the Chandra images, map the spatial distribu- 
tion of the individual temperature components and compare it with 
new Ha images obtained at the Lick Observatory, and archival data 
from the Hubble Space Telescope (HST). In Section 2, we describe 
the data reduction and analysis; in Section 3, we present the results 
of the observations; in Section 4, we determine the physical prop- 
erties of the cooler plasma phases and discuss the uplift of the low 
entropy gas from the bottom of the gravitational potential well, and 
effects of shocks and mixing. Finally, in Section 5, we summarize 
our conclusions. 

A distance to M 87 of 16.1 Mpc ( [Tonry et al.|2001| ) implies a 
linear scale of 4.7 kpc arcmin"^ All errors are quoted at the 68 per 
cent confidence level for one interesting parameter. 
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Figure 1. Left panel: Relative deviations of the surface brightness from a double beta model fit to the Chandra 0.6-2.0 keV image of M 87, obtained by 
coadding all ACIS-I observations (see Paper I for details). The cavities in the core and the two X-ray arms are clearly visible. The XMM-Newton RGS 
extraction region is overplotted. Right panel: The 90 cm radio image by [Owen et al.|{2000) . Two flows (radio arms) emerge from the inner-jet region, one 
directed Eastward (spatially coincident with the X-ray arm) and the other directed to the Southwest. While the Eastern radio arm ends in edge-brightened torus- 
like vortex rings, the Southwestern radio arm develops a S-shaped Southward twist. Both radio arms are immersed in a pair of large, partially overlapping 
radio lobes. Both panels show the same 14.6' x 15.8' field. 



2 DATA REDUCTION AND ANALYSIS 
2.1 Chandra data 

The Chandra observations of M 87 were taken between July 2002 
and November 2005 using the Advanced CCD Imaging Spectrom- 
eter (ACIS). The observations are listed in Table [T] The total net 
exposure time after cleaning is 574 ks. We follow the data reduc- 
tion procedure described in Paper I (see also |Million & Allen|2009[ 
[Million et al.|2009l ). 

The individual regions for the spectral analysis were deter- 
mined using the Contour Binning algorithm ( [Sanders |2006| ), which 
groups neighboring pixels of similar surface brightness until a de- 
sired signal-to-noise threshold is met. In order to have small enough 
regions to resolve substructure and still have enough counts to fit a 
simple multi-temperature model, we adopted a signal-to-noise ratio 
of 50. The total number of regions in this analysis is ~6000, with a 
total number of -15,000,000 net counts. 

Background spectra for the appropriate detector regions were 
extracted from the blank- sky fields available from the Chandra X- 
ray Center. These were normalized by the ratio of the observed 
and blank- sky count rates in the 9.5-12 keV band (the statistical 
uncertainty in the observed 9.5 - 12 keV flux is less than 5 per cent 
in all cases). The background level in these observations is low and 
represents only a marginal source of systematic uncertainty in the 
determined quantities. 

Spectral modeling has been performed with the SPEX pack- 
age ( [Kaastra et al1[T996l SPEX uses an updated version of the 
MEKAL plasma model with respect to XSPEC). We analyze data 
in the 0.6-7.0 keV band. To investigate the multi-temperature struc- 
ture and map the spatial distribution of plasmas at diflferent temper- 



atures, we fit to each bin a model consisting of collisionally ion- 
ized equilibrium plasmas at four fixed temperatures, with variable 
normalizations and common metallicity. This method was intro- 
duced by Sanders et aL[ ( [2004[ ) in the analysis of deep Chandra 
observations of the Perseus cluster. For M 87, the temperatures 
are fixed at 0.5, 1.0, 2.0, and 3.0 keV. The temperatures of the 
coolest and hottest components were selected based on previous 
multi- temperature analysis of XMM-Newton RGS and EPIC data 
( [Werner et al._2006^, Simionescu et al. 2008 ). We also searched for 
a 0.25 keV component in our spectral fits, but did not detect any 
emission at this temperature. The temperatures of the individual 
spectral components are sufficiently far apart that the spectral anal- 
ysis can constrain their emission measures if they are simultane- 
ously present in the same extraction region. The overall metallicity 
in each region is free to vary. The 0/Fe ratio is fixed at 0.59 Solar, 
Ne/Fe at 1.20 Solar, and the Mg/Fe at 0.60 Solar, as determined 
from XMM-Newton RGS spectra ( [Werner et al. 2006 ). The Galac- 
tic absorption toward M 87 is modelled as neutral gas with Solar 
abundances with a column density fixed at A^h = 1.94 x 10^^ cm~^, 
the value determined by the Leiden/Argentine/Bonn (LAB) Sur- 
vey of Galactic Hi (Kalb erla et al.[[2005| ). Throughout the paper, 
abundances are given with respect to the "proto-solar values" by 
[Lodders[p003| ), which are for O, Ne, and Fe approximately 30 per 
cent lower than the values given by (Anders & Grevesse[ ( |1989] ). 

2.2 ZMM-A^^>vron RGS data 

The multi-temperature structure of the hot plasma in and around 
M 87 exhibits a complex set of spectral lines that cannot be mod- 
elled by a single-temperature spectral model ( [Sakelliou et al.[2002{ 
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Figure 2. Spatial distribution of the emission measure, Y = J nuncdV, (in 10^^ cm~^ arcsec"^) of the 0.5 keV (upper left), 1.0 keV (upper right), 2.0 keV 
(bottom left), and 3.0 keV (bottom right) plasma detected at 99.7 per cent confidence. On the upper right panel, we over-plotted the contours of the 90 cm 
radio image jOwen et al.| 2000 >. All the panels show the same 10.45^x9.75' field. For a zoomed in version of the spatial distribution of 0.5 keV plasma in the 
core see the right panel of Fig. [4] 



[Werner et al.||2006] ). In order to confirm that our simple four- 
temperature model describes the temperature structure of M 87 
well, we fit the same model to the XMM-Newton RGS spectra. 

The XMM-Newton RGS data were obtained in January 2005, 
with a net exposure time of 84 ks. The data were processed as de- 
scribed in I Werner et al . | ( [2006] ) . Spectra were extracted from a l.T 
wide extraction region, centred at the core of the galaxy. Because 
the RGS operates without a slit, it collects all photons from within 
the l.Tx ~12' field of view. Line photons originating at angle 
(in arcminutes) along the dispersion direction are shifted in wave- 
length by A/l = 0.1 38A^ A. Therefore, every line is broadened by 
the spatial extent of the source. To account for this spatial broad- 



ening in our spectral model, we produce a predicted line spread 
function (LSF) by convolving the RGS response with the surface 
brightness profile of the galaxy derived from the EPIC/MOSl im- 
age along the dispersion direction. Because the radial profile of a 
particular spectral line can be diflferent from the overall radial sur- 
face brightness profile, the line profile is multiplied by a scale factor 
5, which is the ratio of the observed LSF to the expected LSF. This 
scale factor is a free parameter in the spectral fit. 

We fit the spectra in the 8-38 A band with the same four- 
temperature model fitted to the Chandra data, with an additional 
0.25 keV component included to determine an upper limit on the 
amount of such cool plasma. We have also separately fitted the RGS 
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data with a model consisting of four isobaric cooling flows with 
elemental abundances tied between the models (similar to one of 
the models used for the Centaurus cluster in [Sanders et al.pOOS] ). 
In this case, we model separately gas cooling from 3 keV to 2 keV, 
from 2 keV to 1 keV, from 1 keV to 0.5 keV, and from 0.5 keV 
to 0.25 keV. The central AGN cannot be excluded from the RGS 
data of the core of the galaxy and also needs to be accounted for. 
We model the AGN spectrum with a power-law of a photon index 
7 = 1.95 and a 2-10 keV luminosity of Lx = 1.15 x 10^^ ^^.g ^-i 
( [Werner et allOTg) . 



2.3 Ha data 

Narrow-band Hof imaging (central wavelength 6606 A, FWHM 
75 A) of M 87 was carried out using the Shane 3 m telescope at 
Lick Observatory on March 26th, 2009. The total integration time 
was 24 ks, split into 20 exposures using a dither pattern with ~ 5'' 
offsets. We also acquired 1200 seconds of R band imaging, also 
split into 20 exposures. The data were processed using the pipeline 
of [Erben et al.[ ( [20()5] ). A constant background was subtracted, es- 
timated from the region of the image with the lowest counts. The 
images were registered wit h scamp (Bertin^2006 ), and resampled 
and combined with swarp ( [Bertin et al.|2002| )~The seeing of the 
coadded images is about 2" . 

To study the morphology of the brighter Hof filaments in the 
cluster core, we have also analyzed two Hof images available in 
the HST archive. Both images were taken with WFPC2, placing 
the core onto the Planetary Camera. The Southeastern region was 
imaged with the Wide Field Camera for only 2700 seconds (pro- 
posal ID: 5122, PI Ford), and is not as deep as our ground-based 
image. The Northwestern region was imaged significantly longer 
(13900 seconds, proposal ID: 6296, PI Ford). 



3 RESULTS 

3.1 Spatial distribution of the X-ray gas 

Fig. [2] shows the spatial distribution of the emission measures for 
the four individual temperature components (0.5, 1.0, 2.0, and 
3.0 keV) obtained from the Chandra data. The top left panel shows 
an interesting result: the map reveals the presence of gas with a tem- 
perature of kT ~0.5 keV to the North of the jet in the core of the 
galaxy, in the Southwestern arm, and in a horseshoe like filament to 
the Southeast of the core outside of the mushroom shaped Eastern 
radio arm (radio contours are overplotted in the top right panel of 
Fig.|2]). No 0.5 keV gas with significance higher than 3cr is present 
within the Eastern radio arm. While all the detected 0.5 keV plasma 
in the core and in the Southeastern horseshoe is within the RGS ex- 
traction region (see the left panel of Fig. [TJ, the Southwestern arm 
is outside of the area covered by the present RGS data. 

The spatial distribution of the ~1 keV plasma seen in the top 
right panel of Fig. [2] correlates with the radio emission and is very 
similar to the morphology of the X-ray arms seen clearly in Fig.[T] 
In the core, the 1 keV plasma forms a prominent ridge to the North 
of the AGN, which is cospatial with 0.5 keV plasma and HQf+[Nii] 
emission (see below). Based on two-temperature spectral fits to 
XMM-Newton data it has been shown that the X-ray arms are rela- 
tively isothermal at ~1 keV (Molendi 2 002|[Simionescu et al.|20"08l ) 
and our spectral fits to the Southwestern X-ray arm confirm that its 
temperature is constant as a function of radius (see Fig. [3]). Most of 
the Eastern radio arm contains ~ 1 ke V plasma, a large amount of 




2 3 
Radius (arcmin) 



Figure 3. The temperature profile of the Southwestern X-arm fitted with 
a two-temperature model. While the temperature of the hotter component 
shows a slight radial increase (black data points) the radial distribution of 
the cooler temperature component (red data points) looks remarkably flat at 
around 1 keV. If the uplift were adiabatic, the temperature of the cooler gas 
within the arm would be radially decreasing. 



which is also present to the South of the "stem of the radio mush- 
room" where the 0.5 keV horseshoe is observed. The brightest part 
of the Southwestern radio arm is not spatially coincident with the 
0.5 keV and 1 keV plasma, but it bends around the cooler X-ray 
gas. On the Southwest, just outside of the core region, the 1 keV 
plasma forms a narrow and remarkably straight and smooth fila- 
ment, which broadens with increasing distance from the core. At 
around 3\ which is the approximate distance of the circular shock 
front described in Paper I and by [Forman et al.[ ( [2007| ), the distri- 
bution of 1 keV gas broadens, and at 4.5' it starts to bend to the 
Southeast. 

The spatial distribution of the plasma at ~ 2 keV is, in com- 
parison, remarkably symmetric and only slightly elongated in the 
direction of the X-ray arms. A knot of ~3 keV plasma is present at 
the inner edge of the X-ray cavity South of the core, at the bottom 
of the bright radio cocoon, which Forma n et al.[ ([2005| identified 
as a "bud". This 3 keV plasma is most probably associated with a 
shock front at the radius of 0.6' ( [Forman et al. 12007) Paper I, see 
also the right panel of Fig. [5]). 

The fraction of the emission measure of cooler gas (1 keV and 
0.5 keV) in the X-ray arms with respect to the total emission mea- 
sure is 10-35 per cent, consistent with the results obtained by fitting 
a diflferential emission measure distribution to XMM-Newton data 
by [Simionescu et al.] ( [2008| ). The spatial distribution of the coolest 
gas phases is also in good agreement with that found in the XMM- 
Newton data analysis ( [Simionescu et al.|2 008), but is revealed by 
Chandra with a strikingly better spatial resolution. 

The elemental abundances obtained from the fits are likely to 
be biased due to the limitations of our model, which has its tem- 
peratures fixed at specific values (detailed study of the metallicity 
will be presented in Million et al. in prep.; see also Paper I and 
[Simionescu et al.|200"8| ). 
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Figure 4. The Ha emission (left panel) and plasma at -0.5 keV (right panel; blue contours in the left panel) spatially correlate in the core and form a horseshoe 
like feature to the Southeast of the core. The right panel is a zoomed in version of the top left panel of Fig.|2] Both panels show the same 2.7'xl.85' field. 
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Figure 5. Left panel: Ha image of the innermost Tx2' region of M 87, divided by the best fit De Vaucouleurs profile, obtained by HST. Right panel: Chandra 
pressure map (in units of keV cm~^ x (^MEajl!^^' Paper I for details on producing the pressure map) of the same inner 2'x2' region overplotted with the 
6 cm radio emission contours from Hines et aL| jl98^ . The pressure map clearly shows a discontinuity, a likely shock front, at a radius of 0.6'. The Ha image 
shows that all the bright Ha filaments appear in the downstream region of the shock. The arrows identify the bright Ha features furthest from the centre, at the 
edges of apparently underpressured X-ray cavities, and a bright Ha "knee" on the Southeast. 



3.2 Morphology of the Ha nebulae: new clues about their 
powering 

The left panel of Fig. |4] shows an excerpt of the HQf+[Nii] image 
obtained at the Lick Observatory after subtraction of the scaled R- 
band continuum image. Owing to the large depth, our data reveal 
a previously undiscovered horseshoe-like filament of HQf+[Nii] 
emission to the Southeast of the core. The feature is detected at a 
mean level of 5cr above the sky. Fig.|4]shows that the spatial distri- 
bution of the HQf+[N ii] emission is remarkably similar to that of the 
0.5 keV plasma, both in the core and in the Southeastern horseshoe 
like region. Except for a small knot of Haf+[Nii] emission in the 
"cap of the mushroom" (sec iGavazzi et al.|2000 ), the Eastern radio 
arm is devoid of observable Ha filaments. As mentioned above, the 
Eastern radio arm is also devoid of any detectable 0.5 keV emis- 



sion. Our Lick Observatory data (Fig.|4]) and the HST data (Fig. [5]) 
unfortunately do not cover the Southwestern X-ray arm. 

Fig. |5] shows an Ha image for the innermost Tx2' region of 
M 87 from HST next to a pressure map obtained in Paper I by 
using Chandra data. In the pressure map we see a discontinuity, 
likely shock front, at 0.6' and several underpressured regions filled 
with radio emitting plasma, as shown by the contours of 6 cm radio 
emission from |Hines et al.| ( |1989| ). Assuming that this inner shock 
propagates at an approximately constant velocity of 1000 km s~\ 
it originated in an AGN outburst ~2.8 Myr ago. Interestingly, all 
the bright Ha emission is within the inner higher pressure region 
surrounded by the shock front at 0.6'. Some filaments appear to be 
parallel to the shock front, others seem to lie at the edges of cavi- 
ties, while the filament to the Southeast of the core is perpendicular 



© 0000 RAS, MNRAS 000, 000-000 



AGN heating and gas uplift in M 87 1 



to the shock front, bright in the downstream region and faint in the 
upstream. Upstream of the shock front, the Southeastern Hof fila- 
ment continues at a much lower surface brightness and forms the 
prominent horseshoe. On the Hor image, we identify with arrows 
the bright features furthest from the centre, at the edges of X-ray 
cavities, and a bright Hof "knee" on the Southeast. The same ar- 
rows are plotted in the pressure map. 

3.3 High-resolution XMM-Newton RGS spectra 

To confirm the presence of the coolest gas phases indicated by fit- 
ting the Chandra data and to verify that the simple four-temperature 
model fits the data sufficiently well, we also analyze the XMM- 
Newton RGS high-resolution line spectra. The results of the five- 
temperature fit to the RGS data are listed in Table [2] and the best fit 
model to the spectra is shown in Fig.|6] 

The emission measures for the 2 keV and 1 keV phases ob- 
tained by the RGS are consistent with those from the Chandra maps 
for the area covered by the gratings. The emission measure of the 
0.5 keV phase obtained using the RGS is higher by a factor of 3 
than the value from the Chandra map. The sensitivity of Chandra 
in the 0.6-7 keV band to 0.5 keV emission in these moderate S/N 
spectra is limited and if we only consider regions where the 0.5 
keV component was detected at larger than the 99.7% confidence 
level then we miss the cool emission in some regions. By includ- 
ing all regions detected at the 68% confidence the 0.5 keV features 
in the map become wider and the agreement between the emission 
measures from Chandra and RGS improves. 

The lack of Ovii lines provides strong constraints on the 
amount of 0.25 keV plasma. The 95 per cent upper limit for its 
emission measure is lower by a factor of 11 than the emission mea- 
sure of the 0.5 keV plasma. 

To get constraints on the amount of radiative cooling, we also 
fit the RGS spectra with a set of cooling flow models (see Table [2|. 
If the plasma would cool isobarically in the absence of heating, 
the best fit mass deposition rates would be the same in the diff'er- 
ent temperature ranges. However, the mass deposition rates are de- 
creasing as a function of temperature and for the gas cooling from 
0.5 keV to 0.25 keV we determine a tight 95 per cent upper limit of 
M = 0.06 Mo/yr. 



4 DISCUSSION 

4.1 Orientation of the X-ray arms 

The inferred properties of the X-ray arms depend critically on their 
orientation with respect to our line-of- sight. Based on HST obser- 
vations of the superluminal motion in the jet of M 87, |Biretta et al.| 
( |1999 |) concluded that the position angle of the jet is < 19° from 
our line-of- sight. The two large, partially overlapping, outer radio 
lobes of M 87 (the diff'use structures which embed the more colli- 
mated, brighter radio arms seen in the right panel of Fig. [T] and in 
[Owen et al.|2000| ) are also probably oriented close to our line-of- 
sight. It is likely that the Southern radio lobe is positioned towards 
us and the Northern lobe is oriented away from us. This picture is 
supported by the observation of polarized flux from the Southern 
lobe ( [Andernach et al.|1979| ), which shows that the Faraday depth 
towards it is likely to be small. The radius of the large lobes/bubbles 
is ~22 kpc and their approximate centres are in projection 26 kpc 
apart. Assuming that the physical distance between the centres of 
the bubbles is at least twice their radius, a conservative limit for the 



Table 2. Best fit parameters for a five-temperature fit and a four-cooling- 
flow model fit to the high-resolution RGS spectra extracted from a 1.1' wide 
region centred on the core of M 87. Emission measures, Y = J nunedV, are 
given in 10^^ cm~^. The scale factor s is the ratio of the observed LSF to 
the expected LSF based on the overall radial surface brightness profile. The 
upper limits are quoted at their 95 per cent confidence level. Abundances 
are quoted with respect to the proto-solar values of [Lodders] j2003 ) . 



Parameter 5T-model 4-c.f. model 



Jo.25keV 


< 0.07 


- 


Jo.SkeV 


0.79 ± 0.05 




Mo,5-0,25keV 




< 0.06 


J'l.OkeV 


7.1 ±0.3 




Mi,0-0,5keV 




0.90 ± 0.03 


J^2.0keV 


56.3 ± 1.3 




M2.O-l.OkeV 




3.79 ± 0.25 


Ys.OkeY 


< 2.4 




M3.O-2.OkeV 




6.03 ± 0.21 


S 


2.21 ±0.10 


2.19±0.11 


c 


0.75 ±0.21 


0.82 ± 0.22 


N 


1.9 ±0.3 


2.0 ±0.3 





0.79 ± 0.04 


0.84 ± 0.04 


Ne 


1.92 ±0.13 


1.92 + 0.13 


Mg 


1.25 ±0.15 


1.26 ±0.13 


Fe 


1.33 ±0.04 


1.39 ±0.05 


Ni 


0.9 ± 0.3 


0.9 ± 0.3 




10 20 

wavelength (A) 



Figure 6. The first order RGS spectrum extracted from a l.T wide region 
centred on the core of M 87. The continuous line represents the best fit 
model to the spectrum. 



orientation of the axis along which the bubbles rise is < 35° from 
our line-of- sight. The Southwestern radio arm appears to "merge" 
into the Southern lobe indicating that its orientation is similar to 
that of the two large bubbles. Although the relation of the Eastern 
radio arm to the outer lobe is less clear, the most likely orientation 
of both X-ray and radio arms is approximately anti-parallel and far 
from the plane of the sky. 
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Figure 7. Cumulative gas mass, entropy (S = kT/rr'^), and radiative cool- 
ing time of the 0.5 keV phase (red lines), 1.0 keV phase (blue lines), and 
the ambient ICM (black lines) as a function of projected radius for orien- 
tation angles of 15° (dashed lines), 30° (dotted lines), and 90° (full lines) 
from our line-of-sight. For the deprojected density and temperature profiles 
of the ambient ICM we assume the profiles parametrized by equations (6) 
and (7) in'Churazov et al.|(2008 ). The cooler gas phases are assumed to 
be in pressure equilibrium with the ambient medium (see text for details). 
The vertical dotted lines denote the radii of the outer edge of the ridge of 
multiphase gas North of the jet, the end of the Southeastern horseshoe, and 
the radius of the shock at ~ 160'' which corresponds to the inner edge of 
the "cap of the mushroom" of the Eastern arm. 

4.2 Physical properties and origin of the X-ray arms 

We determine the properties of the X-ray arms for three different 
orientation angles (15°, 30°, and 90°) from our line-of-sight. For 
each orientation, we determine the radial profiles of the gas mass, 
entropy, and cooling time for the different phases (see Fig. [7]). We 
assume that all the gas phases are in equilibrium with the ambi- 
ent pressure at the given distance from the centre. The radial dis- 
tribution of the ambient pressure is determined using the depro- 
jected temperature and electron density profiles parametrized by 
equations (6) and (7) in |Churazov et al.| ( |2008| ). Using this pres- 
sure profile p(r), we calculate the electron number density pro- 
file n^{r) = p{r)/{kT) of the kT = 0.5 keV and 1 keV plas- 
mas. Using these density profiles and the best fit emission mea- 
sures, Y = j nunedV, for each extraction region, we determine the 
volumes which the cooler phases occupy: V = Y/inuno), where 
nn = yiQlX.l is the hydrogen number density for collisionally ion- 
ized plasma with Solar metallicity. 

The line-of-sight depths, /, of the emitting volumes, V, shown 
in Fig. [8] are then calculated as / = V/A, where A is the area on the 
sky of the given extraction region. The corresponding gas masses 
are calculated as M = lAnurripV (where nip is the proton mass), 
only in the spatial regions where their presence was determined 
with better than 3cr significance. The cumulative radial gas mass 
profiles are shown in the bottom panel of Fig. [7] The middle panel 



of this figure shows the radial entropy profiles calculated as S - 
kT/nl^^. The top panel shows the cooling time profiles calculated 
as the gas enthalpy divided by the energy lost per unit volume of 
the plasma: 

|(«e + ni)kT 

^COOl — TT^T^N ' 

f7e^iA(r) 

where the ion number density rii = 0.92;ie and A(T) is the cool- 
ing function for Solar metallicity tabulated by S utherland & Dopitaj 
( 1993|. Cooling functions based on more up to date plasma codes 
(Schure et al.||2009} predict for the 0.5 keV plasma a 9 per cent 
higher cooling rate. 

If the arms lie in the plane of the sky, then the small line- 
of-sight depths of the emitting volumes in Fig. |8] indicate that the 
1 keV gas is not volume filling at most radii, including the core of 
the galaxy. The 1 keV gas is then made of small blobs and filaments 
with an entropy smaller than the lowest entropy of the ambient 
medium in the centre of the galaxy (see Fig.|7]). This would mean, 
as pointed out by Molendi ( 2002 ) , that the cool blobs/filaments can 
hardly originate from the adiabatic evolution of the hot-phase gas 
entrained by buoyant radio bubbles, as suggested by jChurazov etal.| 
(2001 ). However, if the arms are oriented close to our line-of-sight, 
the physical picture changes. 

For the likely angles of 15°-30° from our line-of-sight, both 
the radial distribution of the gas entropy and the distribution of 
the depths of the emitting volumes of the 1 keV component give a 
more consistent picture with other observables. The 1 keV plasma 
appears to be mostly volume filling, except in regions where the 
X-ray images reveal cavities. Its cooling time is >1.5 x 10^ yr, 
which is comparable to or longer than the uplift time. Moreover, 
by taking the projection effects properly into account, the appar- 
ent large angle between the current direction of the jet and the 
Southwestern arm translates into a much smaller physical angle. 
The inferred spatial distribution of the depths of the emitting vol- 
umes (see Fig. [8]) indicates that the orientation angle of the filament 
is not constant and at certain positions the arms are folded along 
our line-of-sight. The Southwestern arm appears to be folded at the 
projected radius of 3.4', where the arm divides into two filaments 
(see Fig.[TJ. The Eastern arm appears to have folds at ~2.r and its 
line-of-sight depth increases in the cap of the mushroom. We cau- 
tion that because of the strongly simplified assumptions the values 
for the depths quoted in the figures are approximate numbers only. 

If the uplift of the low entropy plasma were strictly adiabatic, 
then the temperature of the arms would decrease as a function of 
radius. However, as shown by |Molendi| ( |2002| ), [Simionescu et al.| 
( 2008 ), and in Fig. [3] using a two-temperature model, the average 
temperature of the arms is remarkably constant. The constant tem- 
perature implies rising entropy with radius and suggests that con- 
duction plays an important role in the thermal evolution of the arms. 

The estimated total mass of 1 keV gas in the arms is between 
6-9 X 10^ Mq, which given the uncertainties is remarkably sim- 
ilar to the total gas mass of 6.3 x 10^ M© within the innermost 
3.8 kpc region (0.8' in the plane of the sky), where the X-ray emis- 
sion is relatively spherically symmetric. Outside of this radius, all 
the ~1 keV gas is distributed along the two X-ray arms. If the en- 
tropy of the arms is increasing due to heat conduction, then most 
of the 6-9x10^ M© of ~1 keV plasma might have been uplifted 
from within the innermost approximately spherical 3.8 kpc region 
of the galaxy. While this is a large fraction of the total gas mass 
in the core of the cluster, it is only a fraction of the uplifted gas 
mass observed in more massive clusters; e.g. th e inferred mass of 
the uplifted gas in Hydra A is 1.6-6.1 x 10^ M© ([Simionescu et al. 
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Figure 8. Maps of the line-of-sight depths in kpc for the emitting volumes of 1 keV plasma for orientation angles of 90° (top left panel), 30° (top right panel), 
and 15° (bottom left panel) from our line-of-sight. The map on the bottom right panel shows the line-of-sight depth distribution of the 0.5 keV component for 
an orientation angle of 30° from our line-of-sight. All the panels show the same 10.36^x9.25' field. 



|2009b| ). Though seen with great clarity in our study, the physics of 
gas uplift in M 87 is by no means extreme. Many brightest cluster 
galaxies can be expected to be similarly efficient in evacuating their 
lowest entropy plasma during AGN outbursts. 

By assuming that the outer radio halo is a spherical bubble of 
hot plasma undergoing a steady energy input from the jet, |Owen| 
et al. ( |2000| ) estimated its age to ~10^ yr. However, if the radio halo 



consists of two bubbles rising buoyantly in the hot ICM, then the 
outer radio lobes are l ikely to be older. The a ge of the radio and X- 
ray arms estimated by Churazov et al. ( 200T| is a few times 10^ yr. 



However, if the whole Southwestern arm originated in a single up- 



lift following an AGN outburst, then for an orientation angle of 
<30° from our line-of-sight, its length is >60 kpc, which assuming 
an uplift velocity of 400 km s~^ implies an age of >1.5 x 10^ yr. 
The age of the observed structures is thus greater than previously 
thought. 



4.3 Implications for gas motions in the ICM 

The NW edge of the Southwestern X-ray arm is remarkably smooth 
out to the radius of ~3.3' and the filament is surprisingly straight 
for a distance of over 2\ This strongly suggests that the gas mo- 
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tions cannot be very turbulent in either phase. The lack of strong 
turbulence in the ICM has been inferred previously based on the 
straightness of the Ra filaments in the Perseus cluster jFabian et al.| 
|2003| ), and strong upper limits on turbulent velocities were placed 
based on resonance scattering in NGC 4636 ("Werner et aL"2009 ) 
and turbulent line broadening in Abell 1835 ( Sanders et al. 2010b). 
The bending of the X-ray arms of M 87 at larger radii might be 
caused by gas sloshing, which creates a spiral like shearing motion 
pattern in the ICM which is the most obvious at the two opposite 
and staggered cold fronts seen at 1' to the Southeast ( Sim ionescu| 
let al. 12007 J and at W to the North ( |Simionescu~et al.,2010l ). The 
sloshing induced shearing motions might also generate mixing of 
the ~ 1 keV gas with the ambient ICM at large radii. 



4.4 The coolest X-ray emitting phase 

The total mass of -0.5 keV gas in M 87 is at least 7 x 10^ M©. 
Because we calculate the mass of the 0.5 keV component only in 
spatial regions where it was detected at higher than the 99.7 per 
cent confidence, this value is a conservative lower limit. The in- 
ferred small line-of- sight depths of the emitting volumes of the 
0.5 keV component imply that the coolest X-ray emitting gas is 
not volume filling, but forms a multi-phase medium together with 
both the hotter X-ray emitting and the cold Ka emitting gas phases. 
Soft band X-ray images ( Forman et al.|2007') reveal that this cooler 
phase forms a network of filaments, which spatially coincide with 
Ha filaments ( [Sparks et al.|2004| see also Fig.|4]). Relatively cool 
(~ 0.5 keV) X-ray gas associated with Ka filaments has also been 
observed in the Perseus cluster and in 2A 0335+096 ( [Sanders &| 
[Fabian|2007|[Sanders et an20 09 ). 

One of the most prominent regions in M 87 occupied by multi- 
phase plasma is a ridge to the North of the jet, which contains 2 keV, 
1 keV, 0.5 keV, and Ha gas. This region is just outside of the inner 
radio cocoon and possibly contributes to the confinement of the ra- 
dio plasma on the Northern side and to the depolarization seen there 
by [Owen et al.| ( p^90l ). The depolarization of the radio emission in 
this region overlapping with the filaments indicates that the mag- 
netic fields are disordered on scales less than 0. 1 kpc. Therefore, it 
is likely that the Ha filaments consist of many narrow magnetized 
threads, which are not perfectly aligned and are small compared 
to the observed beam. Depolarization of the radio emission at re- 
gions containing Ha and soft X-ray emission has been previously 
observed in the Centaurus cluster ( ^Taylor et al.|2007| ). 

4.4.1 Limits on cooling and constraints on viable heating 
mechanisms 

Even though the cooling time of the ~0.5 keV gas is relatively 
short (see the top panel of Fig. [7]), no X-ray emitting plasma below 
0.5 keV is observed in M 87. The 95 per cent confidence upper limit 
for the emission measure of 0.25 keV plasma based on the non- 
detection of O VII lines with RGS is 7 = J n^n^dV = 7x10^^ cm~^ 
which is only 9 per cent of the emission measure of the 0.5 keV 
plasma. Therefore, assuming pressure equilibrium between diff'er- 
ent gas phases, the mass of the 0.25 keV phase is at least 22 times 
lower than the mass of the 0.5 keV plasma. Assuming the 1 keV gas 
is isobarically cooling to 0.5 keV, its best fit mass deposition rate 
is M = 0.9 Mq yr~\ which is 15 times larger than the 95 per cent 
confidence upper limit for the mass deposition rate from 0.5 keV to 
0.25 keV. The upper limits on gas cooling below 0.5 keV rely on 
the O VII lines and thus on the O abundance of the cool gas, which 



we assumed to be the same as that of the hotter phase which is 
fairly constant across the core ( Werner et al.|2006| . For significant 
amounts of gas cooling radiatively below 0.5 keV to be present in 
the core, without detectable O vii emission, the O abundance of the 
cooling plasma would have to be unrealistically low. The lack of 
larger amounts of gas radiatively cooling below 0.5 keV implies 
the presence of heating. A similar temperature floor at ~0.5 keV is 
observed in the Perseus cluster, Centaurus cluster, 2 A 0335+096, 
Abell 262, Abell 3581, and HCG 62, ( [Sanders & FabiaK|[2007l 
[Sander s et al. 2008 2009 2010a ). 

The energy emitted per unit volume in a blob of 0.5 keV 
gas, that is in pressure equilibrium with the surrounding 1 keV 
and 2 keV gas, is 4.8 times larger than the energy emitted by the 
neighboring volume of gas at 1 keV and 30 times larger than for 
the gas at 2 keV. Therefore, any heating process preventing the 
coolest gas phase from cooling below the ~0.5 keV temperature 
floor and providing energy uniformally in a volume- averaged way 
would overheat the surrounding gas and is ruled out by the data. 
The mechanism that heats the coolest gas clearly provides more 
energy per unit volume to cooler, denser phases than it does to the 
hotter phases. 

An alternative explanation is that the coolest X-ray plasma 
cools non-radiatively by mixing. [Sanders et al.j ( |2010a[ ) reached 
similar conclusions for a sample of cool core clusters. The coolest 
X-ray plasma in M 87 might cool non-radiatively by mixing with 
the cold optical line emitting filamentary gas, allowing efficient 
conduction between the two phases ( [Fabian et al.|2002[|Soker et aL] 
[2004[ ). As indicated by the fact that all the bright HQf+[Nii] fila- 
ments end at the innermost shock front, this mixing might be pro- 
moted by the shock, which induces enhanced shearing motions, and 
might also power the detected optical and UV line emission (see the 
next section). 

4.5 The nature of the Ha filaments 

The optical emission line filaments in M 87 are likely due to cooled 
gas, enriched by dust from stellar mass loss (S parks et al.|1993] ), 
which has been uplifted by the buoyant relativistic plasma, together 
with the coolest X-ray emitting phase. The total mass of the ^10"^ K 
Ha emitting gas in the filament syst em in M 8 7 is estimated to be of 
the order of 10^-10^ Mq ( Sparks etaL[l993}, which is comparable 
to the mass of ~0.5 keV gas. However, the total mass of cold neu- 
tral and molecular gas based on the upper limits on CO emission 
within the bright Ha emitting regions is only a few times 10^ M© 
( [Salome & Combes|2008[ ). This implies that the Ha emitting gas is 
not a thin skin on the underlying clouds of neutral and molecular 
gas as in the Perseus cluster, wher e the inferred total m ass of the 
cold gas is as high as -4x10^^ M© ( [Salome et al.|2006 i. Recently, 
[Sparks et"aL]p00 9') discovered filaments of C iv emission in M 87, 
with a total power of 1 . 1 x 10^^ erg s~\ arising from gas at a temper- 
ature of ~10^ K, which spatially coincides with the Ha filaments on 
all scales. This detection indicates that the cool and hot gas phases 
are in thermal communication. The emission strengths are consis- 
tent with thermal conduction, which would create this intermediate 
temperature phase at the interface of the hot ICM and the cooler 
gas. However, even though conduction can quantitatively provide 
the energy seen in the form of line emission ( Sparks et al.|2004| , 
it would have to proceed at close to the Spitzer or saturated levels 
to power the filaments. The very thin and long Ha filaments are 
magnetized (see Sect. 4.2 and F abian et al. 2008 ) and the likely ori- 
entation of the magnetic fields parallel to the interface of the cold 
and hot phase would largely suppress the level of conduction. Thus 
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it seems unlikely that thermal conduction alone can explain the op- 
tical line emission. 

All of the bright Hof emission in M 87 is observed in the down- 
stream region of the innermost <3 Myr old shock front, within a 
radius of ~ 0.6'. This suggests that the Hof emission in M 87 is 
somehow related to shocks in the ICM. The passage of a shock ac- 
celerates the intra-cluster medium, but because of the large density 
contrast the relatively cool Hof emitting gas is left behind, result- 
ing in a sudden strong shear around the filaments. It is likely that 
such shear promotes mixing of the cold gas with the ambient hot 
ICM via instabilities. [Ferland et al.| ( [2009 ) recently showed that the 
spectra of optical filaments observed around the central galaxies of 
cooling core clusters can be explained if the filaments are heated 
by ionizing particles, either conducted from the surrounding re- 
gions or produced in situ by MHD waves. By facilitating contact 
between the hot thermal particles and the cold gas, mixing can sup- 
ply the power and the ionizing particles needed to explain the opti- 
cal line emission spectrum. Mixing of X-ray emitting plasma with 
the cold gas also naturally explains the presence of the ~ 10^ K in- 
termediate temperature gas phase. If the coolest X-ray emitting gas 
cools by mixing with the filamentary gas behind the shock front 
and conductively connects with it (e.g. |Begelman & Fabian|1990| 
[Soker et al.|20Q4| l, then its thermal energy will be radiated away in 
the UV band, in the optical HQf+[Nii], and most likely at infrared 
wavelengths in the presence of dust. 

4.6 Implications for galaxy formation models 

The series of recent 'radio' or 'jet' -mode AGN outbursts in M 87, 
revealed in exquisite detail by the Chandra data, by no means rep- 
resent a unique or extreme phenomenon. For example, a complete, 
flux-limited study of eighteen optically and X-ray bright, nearby 
galaxies, including M 87, by |Dunn et aL| ( |2010| ) shows that 17/18 
sources exhibit some form of current radio activity associated with 
the central AGN (see also [Be"st et al.|2005| |Dunn & Fabian|2006| 
[20081 [Sim|2009l ). [Allen al.| ([2006^ also show that the efficiency 
with which the accretion of hot gas powers jets in giant ellipticals, 
including M 87, is remarkably uniform. Our results for the nearest, 
brightest cluster core show that even the current, relatively modest 
level of radio mode AGN activity has been sufficient to remove a 
sizeable fraction of the coolest, lowest entropy gas from the central 
galaxy. This is the material that would have cooled to form stars 
( [Peterson & Fabian|2006[ ). Our results therefore provide direct ev- 
idence that radio mode AGN input has a profound effect on the 
ability of large galaxies to form stars from their surrounding X-ray 
emitting halos. 



• The coolest detected X-ray plasma, with a temperature of 
~0.5 keV, the HQf+[Nii] emitting gas, and the UV emitting gas are 
cospatial, and appear to be arranged in magnetized filaments within 
the hot ambient atmosphere, forming a multi-phase medium. 

• We do not detect X-ray emitting gas with temperatures be- 
low ~0.5 keV. To remain in steady state and at constant pressure, 
the filaments of 0.5 keV gas would require a 5 times higher heat- 
ing rate than the 1 keV gas, and 30 times stronger heating than 
2 keV plasma. Thus, a uniform, volume- averaged heating mecha- 
nism would overheat the surrounding gas and can be ruled out in 
M87. 

• All of the bright Hof and UV filaments are seen in the down- 
stream region of the < 3 Myr old shock front, within a radius of 
~0.6'. This argues that the generation of Ha and UV emission in 
M 87 is related to shocks in the ICM. 

Based on these observations we propose that shocks induce 
shearing around the denser gas filaments, thereby promoting mix- 
ing of the cold gas with the ambient hot ICM via instabilities. This 
process may both provide a mechanism for powering the optical 
line emission and explain the lack of plasma with temperature un- 
der 0.5 keV: 

• By helping to get hot thermal particles into contact with the 
cool, line-emitting gas, mixing can supply the power and the ioniz- 
ing particles needed to explain the optical spectra of the HQf+[Nii] 
nebulae (see [Ferland et al.|20'09l ). 

• Mixing of the coolest X-ray emitting plasma with the cold 
optical line emitting filamentary gas promotes efficient conduc- 
tion between the two phases, allowing the X-ray gas to cool non- 
radiatively, which may explain the lack of gas with temperatures 
under 0.5 keV. 
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5 CONCLUSIONS 

Using a combination of deep (574 ks) Chandra data, XMM-Newton 
high-resolution spectra, and optical HQf+[Nii] images, we have per- 
formed an in-depth study of AGN feedback in M 87. Accounting 
for the fact that the long X-ray and radio "arms" are seen in projec- 
tion with likely angles of ~ 15°-30° from our line-of- sight, we find 
that: 

• The mass of the uplifted low entropy gas in the arms is com- 
parable to the gas mass in the approximately spherically symmetric 
3.8 kpc core, demonstrating that the AGN has a profound effect on 
its immediate surroundings. This result has important implications 
for understanding AGN feedback in galaxy formation models (e.g. 
[Croton et al.|2006l[De Lucia & Blaizot|2007l[Sijacki et al.|2007[ ). 
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